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Abstract 

We consider theoretically a superlattice formed by thin conducting layers separated spatially 
between insulating layers. The dispersion of two coupled phonon-plasmon modes of the system 
is analyzed by using Maxwell's equations, with the influence of retardation included. Both trans- 
mission for the finite plate as well as absorption for the semi-infinite superlattice in the infrared 
are calculated. Reflectance minima are determined by the longitudinal and transverse phonon 
frequencies in the insulating layers and by the density-state singularities of the coupled modes. We 
evaluate also the Raman cross section from the semi-infinite superlattice. 

PACS numbers: 63.22.+m, 73.21. Cd, 78.30.-j, 
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I. INTRODUCTION 



Coupling of collective electron oscillations (plasmons) to optical phonons in polar semi- 
conductors was predicted more than four decade, ago Q, elementally observed using 
Raman spectroscopy in n-doped GaAs |2j and extensively investigated since then (see, e.g., 
jij]). Contrary, the interaction of optical phonons with plasmons in the semiconductor super- 
lattices is much less studied. A two-dimensional electron gas (2DEG) created at the interface 
of two semiconductors has properties which differ drastically from the properties of its three- 
dimensional counterpart. In particular, the plasmon spectrum of the 2DEG is gap less ^| 
owing to the long-range nature of the Coulomb interaction of carriers, u 2 (k) = vpKok/2, 
where vp is the Fermi velocity and k is the inverse static screening length in the 2DEG. 
Coupling of two-dimensional plasmons to optical phonons has been considered in Refs. j^, IsJ 
for a single 2DEG layer. The resulting coupling is usually non-resonant since characteristic 
phonon energies ~ 30 — 50 meV are several times larger than typical plasmon energies. Still, 
hybrid plasmon-optical-phonon modes are of considerable interest in relation to polaronic 
transport phenomena Raman spectroscopy and infrared optical absorption experiments. 

Plasmon excitations in a periodic system of the electron layers have been discussed in 
a number of theoretical papers 0, Q], disregarding the phonon modes. In the present pa- 
per we analyze the coupled phonon-plasmon modes for a superlattice of 2D electron layers 
sandwiched between insulating layers and demonstrate a possibility of stronger resonant 
coupling of plasmons to optical phonons excited in the insulator. This enhancement occurs 
in superlattices due to interaction of plasmons in different layers which spreads plasmon 
spectrum into a mini-band spanning the energies from zero up to the new characteristic 
energy, Vp{kq/ 'd) 1 ^ 2 , where d is the interlayer distance 10]. This value could exceed typi- 
cal phonon frequencies leading to formation of resonant hybrid modes around crossings of 
phonons and band plasmons. 

Usually the coupled phonon-plasmon modes are considered in the so-called electrostatic 
approximation when the retardation effect is ignored and the terms of uo/c in the Maxwell 
equations are neglected in comparison with the terms having values of the wave vector k. 
This is correct if the Raman scattering is studied, when u has a meaning of the frequency 
transfer. Then, u is much less then the incident frequency of the order of 10 1 ~ ck. But if 
we interested in absorption for the infrared region, when uj is the frequency of incident light 
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and corresponds to the optical phonon frequency, u and ck are comparable. In this case, 
the retardation effect must be fully included. 



The plan of the present paper (preliminary results were published in Ref. ll|) is the 
following. In Sec. II, the Maxwell equations for the periodic system of thin conducting 
layers sandwiched between the insulating layers are solved to find a spectrum of coupled 
phonon-plasmon modes. In Sec. Ill, we consider absorption of the finite sample of the 
layers and reflectance of the semi-infinite system in the infrared region. In Sec. IV, we 
analyze the Raman light scattering from the semi-infinite system of layers. 



II. SPECTRUM OF COUPLED MODES 



We consider a superlattice formed by periodically grown layers of two polar semicon- 
ductors ( e.g., GaAs and AlGaAs) with 2DEG layers formed in the interface regions, see 

Fig.m 



Nd+z 



FIG. 1: The stack of 2d conducting layers separated between dielectric layers of thickness d with 
dielectric constant e(w); zq and Nd + z\ are the boundaries (dashed lines) of the sample if the finite 
stack is considered. 



For the sake of simplicity we assume a superlattice with a single period d and the thickness 
of a 2DEG layer much less than the period. We also neglect the difference in bulk phonon 
properties of the two materials. Optical phonons in polar semiconductors are most conve- 
niently described within the dielectric continuum model, which yields the familiar phonon 
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contribution into the dielectric function, 



e{u>) 



CJ^O — UJ Z — 1001 



where co>lo an d ^to are the frequencies of longitudinal and transverse optical phonons, 
respectively, and V is the the phonon width (we do not distinguish the widths of the TO 
and LO modes). 

Collective modes of our system are described by the Maxwell equations, which in the 
Fourier representation with respect to time have the form 

V(V.E)-V 2 E = 4 E+ ^j, (!) 

where the last term takes into account the in-plane electric currents j induced in the 2DEG 
layers by the electric field E. As usual, when the frequency of collective mode lies above 
the electron-hole continuum, it is sufficient to use the Drude conductivity to describe the 
in-plane electric currents, 

■ o 

j\\(u;,x,z) = — — r^rl^O ~ ^n)E||0,£,z), (2) 
m (u + vy) n 

where z n = nd are positions of interfaces (n is integer corresponding to the periodicity of 
the stack), n e = p 2 F /2-Kh 2 and m are the electron density (per the surface unit) and the 
effective mass, 7 is the electron collision frequency, x and z are the coordinates along and 
perpendicular to the interfaces, respectively. 

We consider the case of the p— polarization, when the field E lies in the xz— plane and, 
therefore, the current j has then only the x— component. Making use of the Fourier transfor- 
mations with respect to the x— coordinate, E ~ e lkxX , we can rewrite the Maxwell equations 
(Q) in the form 

, dE z d 2 E x , Amiuj . 

ik x —i TV - e(u c) E x = — —j x 

dz dz z cr 

dE 

ik x ^ + {k 2 x - euo 2 /c 2 )E z = 0. 



Eliminating 

we get for E x the equation 

[dz 2 



k 2 dz ' 



k 2 - 2kCJ2$( z - z n) ) E x (u, k x ,z) = 0, (3) 

n ) 



where 

2im e e 2 K r— — — 

6 = — t — — , k = \k 2 - ecu 2 c 2 . 

At the interfaces z = z n , the E x component must be continuous and the z— component of 
electric induction eE z has a jump: 

rnd+O 

e(E z \ z=nd+ o - E z \ z=nd - ) = 4vr / p(u, k x , z)dz, (4) 

Jnd-0 

where the carrier density is connected to the current © according to the equation of conti- 
nuity 

p(w, k x , z) = j x (u, k x , z)k x /u. (5) 

For the infinite stack of layers, — oo < n < oo, independent solutions to Eqs. © - © 
represent two Bloch states E x (z) = f±(z) moving along positive and negative directions of 
the z-axis, 

f ± (z) = e ±ifcznd sinh k{z - nd) - e Tik * d sinh k[z - (n + l)d], nd < z < (n + l)d (6) 

with the quasi-momentum k z determined from the dispersion equation: 

cos k z d = cosh nd — C sinh nd. (7) 

The quasi-momentum k z can be restricted to the Brillouin half-zone < k z < ir/d, if 
the phonon and electron damping vanishes. In the general case, we fix the choice of the 
eigenfunctions in Eq. © by the condition lmk z > 0, so that the solution / + decreases in the 
positive direction z. The equation © implicitly determines the spectrum u±(k x , k z ) of two 
coupled plasmon-optical-phonon modes shown in Figs. |2Eas functions of the in-plane wave 
vector k x and the quasi-momentum k z . These modes are undamped if the electron collision 
rate and the phonon width are small. The modes arise from the interaction of the plasmon 
branch in the 2DEG and the phonon LO mode in the 3-d insulator. They have a definite 
character far from the intersection of the corresponding dispersion curves. For instance, at 
small values of k x , the a;_(fc x , k z ) mode is mainly the plasmon mode, whereas the u + (k x , k z ) 
mode has mainly the plasmon character. At large values of k x , they interchange their 
character. 

Note, that k x = k z = is a saddle point for both branches u±(k x , k z ). In the vicinity of 
this point, the frequency grows as a function of k x and decreases with increasing k z . This 
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1.57 



quasi-momentum 



FIG. 2: The plasmon-like (at small wave vectors) mode of an infinite number of 2D metallic layers 
(with the carrier concentration in every layer n e = 3 ■ 10 11 cm~ 2 ) sandwiched between dielectric 
layers of thickness d with the phonon frequency wlo The frequency uj-(k x ,k z ) in units of wlo 
is plotted as a function of the in-plane wave vector k x and the quasi-momentum k z < n/d (both 
in units 1/d). Parameters used here are reported in the literature for GaAs: cjlo = 36.5 meV, 
wro = 33.6 meV, £oo=10.6; the thickness d is taken as d = 1/kq, kq = 2me 2 /?i 2 e 0O = 2.5 • 10 6 
cm -1 is the screening length. 



FIG. 3: Same as in Fig. EJbut for the phonon-like (at small wave vectors) mode u+(k x , k z ). 

is evident in Figs. EEand can be explicitly shown in the limit of k x S> u\e(u)\/c, when the 
retardation effects of the electromagnetic field are negligible. Then, k ~ k x and equation 
(JJJ) yields two solutions, 




quasi-momentum 



ool(k x , k z ) = \{tf + ul ) ± l - [(Q 2 + ut f - 



1/2 



(8) 



where we introduced the notations, 



Vt (k x , kz) — (jJ()(k x ) 



sinh k r d 



2im e e 2 k, 



•x 



cosh k x d — cos k z d 



me. 



oo 
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The frequency uo(k x ) is the conventional square-root spectrum of two-dimensional plasmons 
in the limit of layer separation large compared to the wavelength, d 3> k~ l . The quantity 
Q(k x , k x ) describes the plasmon spectrum that would exist in a superlattice of consisting 
non-polar semiconductors (when uj^q = c^to)- 
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FIG. 4: Dispersion of coupled phonon-plasmon modes of an infinite number of 2D metallic layers 
sandwiched between dielectric layers. The frequencies uj±(k x ,k z ) (in units of wlo) are plotted as 
functions of the in-plane wave vector k x for two values of the quasi-momentum k z (all in units of 
the inverse period 1/d = 2.5 • 10 6 cm -1 ). The dashed-dotted lines represent the upper boundary 
iv±(k x ,k z = 0) of two phonon-plasmon modes, while the solid lines mark the lower boundary 
corresponding to uj±(k x , k z = ir/d). Parameters are the same as in Figs. El -El 
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FIG. 5: Dispersion of coupled phonon-plasmon modes. The frequencies u}±(k x , kg) are plotted as 
a function of the quasi-momentum k z for two values of the in-plane wave vector k x . 

The mode u> + (k x , k z ) has a gap while the other mode has a linear dispersion u-(k x , k z ) = 
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s(k z )\k x \ at k x <C (d 1 , fc*) with the velocity given by 



K()d 



2(1 — cos fe^cf) ' 

where kq = 2me 2 /h 2 e oc is the static screening radius in the 2DEG. For a fixed value of the 
in-plane wave vector k x , every mode develops a band (with respect to the quasi-momentum 
k z ) with the boundaries, oj upper {k x ) < u(k x , k z ) < u)i ower {k x ) where 

k d 

^Ipper = V 2 (k x , 0) = U%(k x ) COth -y- 

for the upper boundary and 

k d 

Slower = M 2 (k x , ir/d) = u 2 (k x ) tanh -y- 

for the lower boundary Fig. El illustrates the behavior of the gapless, so called "acoustic" 
mode for the small values of k z < k x , where it acquires a gap. We see that the frequency of 
this mode decreases rapidly in the region k z > k x . 

In the rest of the paper we analyze various experimental implications resulting from the 
existence of hybrid plasmon-phonon modes. Such modes can be observed in both the infrared 
absorption and the Raman spectroscopy. 



III. INFRARED ABSORPTION ON COUPLED PLASMON PHONON MODES 

We calculate now the reflectance and the transmission of a plane wave with the p- 
polarization, incident from the vacuum on the thin plate consisting of the system of layers. 
Let us suppose that the boundaries of sample are parallel to the layers and intersect the 
z— axis at z = z and z = Nd + Zi, respectively with < z , Z\ < d (see Fig. |TJ). We assume 
in the vacuum (z < z ), 

E x (z) = e ifc ' l) ( 2 - 2 o) + Ae- ik * )(z - Zo \ 

where = J (<^/c) 2 — k 2 and A is the amplitude of the reflected wave. In the region 
z > Nd + z\ , the transmitted wave has the form 

E x {z) =Te tk * ){z - Zl - Nd \ 

and we search inside the plate the field as a sum of two solutions ©: 

E x (z) = C + f+(z)+C-f-{z), 
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FIG. 6: - Fig. Calculated p-polarized reflection-absorption spectra of GaAs plates (of the 
thickness I = 300<i) with the super-lattices of different electron concentrations in a layer; the 
frequencies, the phonon width T, and the electron relaxation frequency 7 are given in units of lulo- 
The lattice period d = 1/ko and the incidence angle 9 = 7r/4. Parameters used here are reported 
in the literature for GaAs: wlo = 36.5 meV, wto = 33.6 meV, £00 = 10.6. 
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with the definite values of k x , co, and k z = k z (u, n) obeying the dispersion equation Q. 

As usual, the boundary conditions at z = zq and z = Nd + z\ require the continuance 
of the a;— component of the electric field E x parallel to the layers and of the z— component 
of the electric induction D z normal to the layers. The second condition rewritten via the 
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FIG. 8: 
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FIG. 9: 



electric field E x gives, e.g., at z = z 

e 1 
—E' x (z +) = -—^(zq-). 

K kl 

The boundary conditions give for C + , C_, A, and T the following equations 



l + A = C + f + (z ) + C^f4z ), 

-l + A=[c + f' + (z ) + C.f'_(z )_ 
T = C + f + (z 1 ) + C.f.(z 1 ), 
-T = 



IK* 



C+f+(z x ) + C-f_(z x ) 
10 



Solving these equations, we find the transmitted amplitude 



T 



2efc«/;(zi)/_(z 1 )-/ + (z 1 )/i(zi) 



in 2 ana 2 2 — 012021 



and the reflected amplitude 



.4 



-1 + 2 



tW-Oo) - «22,/+(>o) 



0-110-22 — 0-120-21 



where 



a u = f' + (z Q )ekV/m 2 - f + (z ),a 22 = f_{z x )ekf /in 2 + /_(z x ), 
012 = t{zo)ekf/i^ - /_(zo),a 2 i = f+(z x )eh® /iJ + / + (z x ). 



The transmission |T| 2 and reflection |A| 2 coefficients are shown in Figs. EHfor samples of 
various electron concentrations n e . The incidence angle is taken 9 = 71 /4. Other parameters 
are the following: the thickness I = 300 d, the lattice period d = 1/kq = 4 • 10 -7 cm, 
the phonon width F = 0.02u;lo, and the electron scattering rate 7 = O.OIculo- in all that 
case k z d << 1, and there are many layers on the wave length of the field. Therefore, the 
reflectance is not sensitive to the positions of the sample surface z and Z\. But, if the 
thickness d is more larger, as shown in Fig. El an interference phenomenon is seen. 



FIG. 10: Interference in the sample of thickness I = 3000<i; the electron concentration is 10 
cm . 

In order to avoid the interference effect, we calculate the reflectance for the semi- infinite 
sample. The results can be seen in Figs. El- [3 where the theoretical curves are presented 
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FIG. 11: Calculated p-polarized reflection-absorption spectra for the semi- infinite super lattice via 
the frequency in units of wlo a t the incident angle 7r/4. The electron concentration is labelled on 
curves, the lattice period is d = 1/kq = 4 x 10 -7 cm; the phonon width is T and the electron 
relaxation rate is 7 (in uns. colo)- 
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FIG. 12: Same as in Fig. II 1| but for larger period. 

for the various electron concentrations and the lattice period. There is a singularity at 
cg>lo F° r d — 1/kq and the intermediate electron concentration in the layer (dashed line in 
Figs. ITTlll2j) . there are two regions (one is at higher frequencies and another is just below 
ujto= 0.9c<jlo) where the sample with layers is more transparent than the sample without 
any (shown in the dashed-dotted line). 

This is an effect of the coupled phonon-plasmon modes: the minima of the reflection 
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coefficient correspond to the density state singularities of the phonon-plasmon coupled modes 
at k z = and k x — > (see Fig. HJ). The frequencies of these singularities are determined by 
Eq. (jSJ) with Q 2 = 47m e e 2 /mGfe 0O . For the large electron concentration (upper solid line), 
the reflection is not complete only in the narrow interval bounded by the singularities at uj^o 
and c^lo Finally, the reflection coefficient tends to the unity at the low frequencies because 
the skin depth of the metallic system goes to the infinity in this case. For the sample with 
the large period d = 5/k,q (Fig. IT2j) . the effect of carriers is seen at the higher concentration. 

IV. RAMAN SCATTERING FROM COUPLED MODES 

Consider now the Raman scattering of the radiation incident from a vacuum with the 
vector potential A*, the frequency u/, and the wave vector k l on the sample occupying the 
semi-infinite space z > z , where < z$ < d. The corresponding quantities in the scattered 
wave are designated as A s , u s , and k s . 

In addition to these fields, an electric field E is excited in the Raman light scattering in 
polar crystals along with the longitudinal optical vibrations u. The field E corresponds to 
the excitation of plasmons, whereas the vibrations u associate with the phonon excitations. 

These processes can be described with the help of the effective Hamiltonian 

H = J d 3 rAf jk (t, v)A]{t, v)A\{t, r), (9) 

where the operator 

N jh (t, r) = g^ k u t (t, r) + g? jk Ei(t, r) (10) 

is linear in the phonon u and photon operators E. 

We are interested in the inelastic scattering on the phonon-plasmon coupled modes. 
Therefore, we assume that the frequency transfer u = uo l — uo s is of the order of the phonon 
frequencies co>LO) but the frequencies u l and uj s of the incident and scattered fields are much 
greater than the phonon frequencies. Then, we can ignore the effect of carriers in the layers 
on the propagation of both scattered A"(t,r) as well as incident A l (t,r) lights, taking the 
incident field in the sample (z > z ) in the form 

4(a/,4,z) = t je ikl{z - zo) , 

normalized to the incident flow; tj is the transmission coefficient from vacuum into the 
sample, k l z = e{ui' l ){ui i / c) 2 — k t2 . 
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When the scattered field A s is taken as the variable in the Hamiltonian Q, we obtain 
an additional current j s (t, r) 

~ c JXs= jj(t> r ) = -oSTjkit, r)A{(t, r) (11) 

in the Maxwell equation for this field. The phonon and plasmon fields in J\fjk(t,r) are 
the source of the Raman scattering, whereas the incident field A*(t, r) is considered as the 
external force for scattering. 

Eliminating the z— component, we write the Maxwell equation for the p— polarization of 
the scattered field (x — z is the scattering plane) in the form 

( d 2 ^ 



(j 2 +k?(cu s )\ E s x (u s ,K,z) = I(u s ,K,z), 



/(- s , K, z) = (a*V)£K K, *) + iK dj ' {w '^ z) ) (12) 



where 



and k s z (uj s ) = ye(uj s ) (lu s / c) 2 — k s x is the normal component of the wave vector in the medium 
for the scattered wave. 

The scattered field in the vacuum z < z is expressed in terms of I(uj s , k x , z): 

xK x ' k s z + e(uj s )uj s cos6 s /c 

E s z (u s ,k xl z) = {k s u a cose 3 /ck?)E s x {u s ,k s x ,z), 

where 

j roc 

h = — dz'e ik ° iz -o)coseVcj^ S) k s^ ^ ( 13) 

and 9 s is the propagation angle of scattered wave in the vacuum. 

The energy flow from the surface of the sample is given by \E x (u s , k x , z < zq)\ 2 . Therefore, 
we have to calculate 

<P(Lu s ,k s x ,z)I(u s ,k s x ,z')> (14) 

averaged quantum-mechanically and statistically, where according to Eqs. ()12j) and (fTT|) we 
meet the Fourier transform of the correlation function 

K jkij , y (t, r; t', r') =< M*(t, r)Af fk ,(t', r') > . 
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Because this correlator depends on the differences (t — t') and (x — x'), we can expand that 
in the Fourier integral with respect to these differences. Then, we have the Fourier trans- 
form Kjk .j'k'{^, k x , z, z') which can be expressed in terms of the generalized susceptibility 
according to the fluctuation-dissipation theorem: 

2 

K jkJ < k ,(uJ,k x ,z,z') = - _ exp( _^ T ) lm Xjkj'k'{^, k x , z, z'). 
The generalized susceptibility Xjk,j'k'{w, k x , z, z') is involved in the response 

< Njk(u:, k x , z) >= - / dz'xjkj>k'(u, K, z, z')Uj>k>{u, k x , z') (15) 



■-o 



to the force 



U jk (u =J- uo% k x = fc* - k s x , z) = A S j(uj s , k s x , z)A\(uj\ 4, z) 



where q z = k\ + k s z . 

To calculate the generalized susceptibility Xjk,j'k'{u, k x , z, z'), we write the equations for 
the averaged phonon u and plasmon E fields. The variation of the Hamiltonian Q with 
respect to the vibrations u gives in the right-hand side of the motion equation for the phonon 
field 

(u4o - ^ 2 - iujT)ui(uj, k x , z) = —Ei(u, k x , z) % -^-U jk (uj, k x , z), (16) 

the additional term to the force from the electric field; p is the density of the reduced mass 
and Z is the effective charge. 

The equation for the plasmon field E can be obtained if this field is taken as the variable 
in the Hamiltonian 0: 

V • (E + 4vrP) = 4vrp. (17) 

The charge density p is connected to the current by Eq. (JHJ). The polarization P includes the 
dipole moment Zu, the contribution of the filled electronic state XooE, and the variational 
term -5H/5E: 

Pi = Zui + XooEi- gf jk U jk . 

Here we can substitute using Eq. (fTtjj). We obtain in Eq. (JTTJ) the term with = l+4nXoo 
and dfjk 9fjk- 

~9%k = 9m + g^z/pi^o ^r). (is) 
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Equation (fT7|) can be substantially simplified in the case under consideration when the 
the wave vectors k*, k s , and consequently the momentum transfer k = k* — k s are determined 
by the frequency of the incident optical radiation, whereas the frequency transfer u of the 
excited fields E and u is much less than the incident frequency, uj = u l — uj s <C uj 1 . Therefore, 
we can neglect here the retardation terms u/c compared to k x (so that k = k x ) and introduce 
the potential, E = — V0. We obtain for the potential the following equation 

~ k 2 x - 2k x Cj2S(z - zM <f>(u, k x , z) = J^j (ik x ~g% k + 9% k ^j U jk (u, k x , z). (19) 

The solution to this equation is found by using the Green's function obeying the same 
Eq. (|T9"j) but with the 5(z — z') function on the right-hand side. As very well known, the 
Green's function is expressed in terms of the solutions of the corresponding homogeneous 
equation. Taking into account f±(uu,k x ,z) (JHJ), we write: 

G(2 , ;<) = ' — ( > ' (20, 

2/«sin (k z d) smh (red) [ f_( z )f + ( z >), z < z > 

where where the quasi-momentum k z = k z (uj s , k) has to be determined from the dispersion 
relation (JJJ) ignoring the retardation, k — k x . 

According to T d symmetry of GaAs lattice, the Raman tensors g^ k and gfj k have only two 
independent components g xxx and g xyz in the crystal axes. Let the scattered light is always 
polarized in the xz— plane. Now we consider two geometry. For the parallel scattering 
geometry (a), the incident light considered to be polarized in the x— direction, then the 
a;— components of exited phonon and plasmon fields can be active in the Raman scattering 
due to g xxx . For the crossed geometry (b), the incident light is polarized in the y— direction. 
Then, z— components of the exited fields are active and we take into account terms with g xyz . 
Thus, for the generalized susceptibility Xxx,xx(u, k x , z, z') defined by Eq. (fTKj) . we obtain 



n u2 In dirk 2 

Xxx>x M kx, z, z') = 2 9xx i P . r 5(z - z') - —^~g° xx G(z, z>) (21) 

Ct/rpQ — UJ — IUJL £[Lu J 

To get the Raman cross-section, Eq. ((HJ), we evaluate the integral 

/•OO 

/ dzdz'e^ z ^^Imx^^.w.z^'). 
We obtain for terms with G(z, z') 

Int= f°°dzf-(z) l°° dz'f + {z')e i{q * z '- q * z) + f° dzf+(z) f ' dz' ' f_{z')e i ^ z '- q >\ 



2Q J Z J 2Q 
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and 



dzf-(z) / dz' f + (z')e l{qzZ ~ q * z) 



00 rz +d rz+d 

e i{q*-q%)nd+i{k z -q z )md dzf_{z) dz' f + (V)^ 2 ~ q * z) . (22) 

n,m=0 J z 

The sum over the integer n could be extended to the infinity being equal to 5/ 2d, if the 
thickness of the sample is larger compared with the skin-depth 5 of the incident or scattered 
waves. The sum over m reduces to a not vanishing factor only under the Bragg condition, 
k z — q z = 27m/ d, which expresses the momentum conservation law in the exciting processes 
of the phonon-plasmon coupled modes. In the macroscopic limit, when the wave length of 
the exited mode is large compared to the period d, only the main Bragg maximum (n = 0) 
is observed for each of the coupled modes. In this case, we can expand in powers of k x d and 
k z d the matrix element in Eq. (|22|). 

Omitting the common factors, we have the Raman intensity for the parallel geometry (a) 
in the form 

Int xx (u, k x ) = Im j 2 9 **{ P . - + (23) 
( e , 9 xxx Z/p V 4vrA;2 ] 

I y xxx ' 



and for the crossed geometry (b) 



Int xy {u, k x ) = Im { 2 9 * 2y { P . - + (24) 
9x yz Z/p \ 2 4vrg 2 2 1 



E I Jxyz 



9xyz ~^ 



The wave- vector k z of the coupled phonon-plasmon mode is determined by Eq. ((7j). For 
example, if the incidence is normal to the sample surface and 9 is the scattering angle, then 



k x = LO l sax9jc and q z = (ye(uj l ) + \J e{uj % ) — sin 6)u l /c, where we take into account that 
uj 1 ~ uj s . 

In Eqs. (12 3p - (|24)1. we dropped slowly varying factors depending on the parameters of 
the incident and scattered radiations, for instance, their dielectric constants as well as the 
temperature factor 1/[1 + exp (—u/T)]. In numerical calculations, we used the relationship 
between the vertices g u and g E known from experiments and given by the Faust-Henry 
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FIG. 13: Raman intensity vs. the frequency transfer (in units wlo) i n the z(xx)z' geometry for 
two values of the current concentration indicated at the curves in units of 10 12 cm -2 ; z and z' are 
the propagation directions of the incident and scattered lights, respectively; the direction z' makes 
an angle of 7r/4 with the backscattering direction — z, .(xx). are the corresponding polarizations. 
The values of the other parameters are the same as in Fig. |1J 
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FIG. 14: Raman intensity vs. the frequency transfer (in units lulo) m the z(xy)z' geometry for 
the angle of scattering varying from (bottom) n/2 (top) with a step of n/10. 

constant Cfh = 9 U ^ / g E pui\ = —0.5. The results of calculations are shown in Figs. Upland 

m 

Note, that for the case of normal propagation of both incident and scattered radiation, 
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k x = 0, the second term in Eq. (|23|1 vanishes, and the Raman peak at parallel polarizations 
(geometry a) is situated at the frequency cuto The others peaks in Fig. H21 correspond to the 
excitations of the coupled phonon-plasmons. But for the crossed polarizations (geometry b) 
and k x —>■ 0, the dispersion equation (J7J) with C —>■ gives k z = ik x . Then, using Eq. (J24j) 
and the relationship u^q — = 4iiZ 2 /e 00 p between the frequencies of the longitudinal and 
transverse phonons, we see that the peak appears only at u;lo, since the terms with poles at 
ujto are cancelled. This peak corresponds to the zero of dielectric constant s{uj). At other 
scattering angles, a peak appears at cjto independently of the scattering angle. Two other 
peaks on each curves in Fig. correspond to the excitation of the phonon-plasmons. 

V. CONCLUSIONS 

In this work, we investigated the infrared absorbtion and the Raman scattering on the 
coupled phonon-plasmon modes with the help of a simple model of superlattices formed by 
thin conducting layers separated with insulating layers. This model admits the dispersion 
relation of an analytical form. Our results for the reflectance and the Raman spectra show 
that the observed picture can be drastically modified by means of the carrier concentration, 
the superlattice period, and the frequency. 
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